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Abstract 

We have evaluated the angular-momentum distribution functions for finite 
numbers of electrons in Laughlin states. For very small numbers of elec- 
trons the angular-momentum state occupation numbers have been evaluated 
exactly while for larger numbers of electrons they have been obtained from 
Monte-Carlo estimates of the one-particle density matrix. An exact relation- 
ship, valid for any number of electrons, has been derived for the ratio of the 
occupation numbers of the two outermost orbitals of the Laughlin droplet and 
is used to test the accuracy of the Monte-Carlo calculations. We compare the 
occupation numbers near the outer edges of the droplets with predictions 
based on the chiral Luttinger liquid picture of Laughlin state edges and dis- 
cuss the surprisingly large oscillations in occupation numbers which occur for 

angular momenta far from the edge. 
PACS numbers: 73.20Dx,73.dO.Mf 



1 



Typeset using REVTgX 



In the strong magnetic field limit all electrons in a two dimensional electron gas occupy 
only the single-particle orbitals which have the minimum quantized kinetic energy, those 
of the lowest Landau level. In the symmetric gauge, appropriate for a finite droplet of 
electrons with circular symmetry, the lowest Landau level orbitals [jXJ] are labeled by angular 
momentum: <f) m (z) = (27rm!2 m ) -1 / 2 z m exp(— |z| 2 /4), m = 0,1,2,... where z% = Xi + k/j. 
(We use the magnetic length, i = (hc/eB) 1 ^ 2 , as the unit of length.) The orbital with 



angular momentum m has its radial coordinate localized within ~ £ of R m = \f2mt while 

its dependence on its angular coordinate is the same as that of an orbital with momentum 
k m = 27Tm/R m in a one-dimensional (ID) electron gas with periodic boundary conditions. 
The system is therefore strongly analogous to an interacting ID electron gas except that 
the orbitals have only one sign of momentum and the effective 'length' of the ID gas system 
depends on the momentum of the orbital. For orbitals near the edge of a large droplet R m 
is nearly fixed and the analogy with a one dimensional gas with only one sign of momentum 
(a chiral Luttinger liquid) can 0|| be made precise. 

In the ground state of an iV-particle non-interacting electron gas the orbitals at energies 
below the Fermi energy are occupied and all others are empty. For an interacting gas the 
single-particle state occupation numbers fluctuate. The change in average momentum state 
occupation numbers, the momentum distribution function, is one of the qualitative changes 
due to electron-electron interactions. In this Rapid Communication we report on the first 
study of the angular-momentum distribution function (AMDF) for strongly correlated two- 
dimensional electron droplets in a strong magnetic field 0. The most compact iV-electron 
droplet is a single Slater determinant in which the orbitals with m = 0, 1, ■ ■ ■ , N — 1 are 
occupied; the area of this maximum density droplet is A^ DD ~ N2tc£ 2 ~ ttR%. When 
the electron system occupies a larger area the average occupation numbers will be smaller 
than one and the ground state of the electron system will be strongly correlated. The best 
understood strongly correlated states are those discovered by Laughlin 0: 





(1) 
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is the ground state when the droplet occupies @-f§ an area ~ MAfJ DD . The M=l Laughlin 
state is the maximum density droplet. We restrict our attention here to these many-body 
states. 

We evaluate the one-body density matrix, n(z, z') defined by 

n(z, z') — N J d 2 z 2 ... J d 2 z N ip M (z, z 2 , ...z N ) x 

tfj* M (z' J z 2 ,...z N )/Q N (2) 

where Qn is the normalization integral for the iV-electron Laughlin wave function. The 
second quantized form of Eq. (0) relates n(z, z') to the AMDF: 

M(N-1) 

n(z,z')= J2 ( n rn)<P* m (z')(i) m (z), (3) 
m=0 

To obtain Eq. (^) we have noted that no orbital angular momentum larger than M(N — 1) is 
ever occupied in the iV-electron Laughlin droplet and that the Laughlin state is an eigenstate 
of total angular momentum. (Mtot = MN(N — l)/2). Eq. (|3|) can be inverted to express 
(n m ) in terms of n(z, z'): 

(n m ) = j d 2 z J d 2 z'(j)* m (z)n(z,z')(i) m (z'). (4) 

For sufficiently small iV and M the AMDF for a Laughlin droplet can be evaluated analyt- 
ically. Some results for very small droplets are listed in Table |. 

For larger droplets we have been unable to evaluate the AMDF analytically [[Kj but 
some properties of the distribution function are known. One such property follows from the 
expansion of iPm[z) in decreasing powers of Z\. 

'(zi, • • • , z n) = 

[ Zl M(N-l) _ M{N _ l)z M(N-l)-l 2 + . . .] x 

exp (-\z 1 \ 2 /4)ilj < y~ 1) (z 2} • • • , z N ). (5) 

In Eq. (^|), ip^~ 1 \z2 } ■ ■ ■ , zn) is the N — l electron droplet, Z is its center of mass coordinate 
and only the two highest powers of Z\ have been retained. Using Eq. (J|) in Eq. (f|) we obtain 



(n mo ) = iV%^2 m °m ! (6) 

where m = M(N—1) is the angular momentum of the outermost occupied orbital. Similarly 

(n mo _x) = iV%^2 m °- 1 (m - 1)\M 2 (N - lf{Z 2 ) N ^. (7) 

It is easy to show a Laughlin droplet has definite center-of-mass angular momentum equal 
to zero || from which it follows that (Z 2 )at-i = 2£ 2 /(N — 1) and hence that for any iV 

(n mo _i) = M{n mo ). (8) 

This result can be extended to orbitals farther from the outer edge using an argument 
due to Wen 0. The density, n(z) = n(z,z) far outside the droplet may be determined 
up to a constant from Eq. (§) and Eq. (||) by using a plasma analogy f5|. The result for 
r = \z\ > i? mo is [|] 

n(r) oc exp(-r 2 /2)(r 2 /2) m °(l - R 2 m Jr 2 )- m . (9) 

Expanding the right hand side of Eq. ( |^) and comparing with Eq. ( ^) gives 

. , (k + M-l\. , 

{n mo -k} = I M _ 1 I {n mo ) (10) 

which reduces to Eq. (§) for k = 1. We remark that while Eq. @ is exact for any number of 
particles Eq. (|10D (for k > 1) becomes exact only in the limit N 1 ^ 2 ^> k. (For example, we 
see that the occupation numbers listed in Table | do not satisfy Eq. (0). Eq. (|T0D implies 
that for N 1 / 2 > k > M, ( n m -k) oc k M 1 . The behavior of the AMDF in this limit is 
thus consistent with Luttinger liquid |L2] behavior with a critical exponent related to the 



quantized Hall conductance as argued on more general grounds by Wen [f2[. This property 
of the AMDF presumably holds as long as the fractional quantum Hall effect occurs and is 
not unique to the Laughlin state. 

To examine the AMDF in the interior of the droplet and to test how well Eq. ( ]TTJ| ) is 
satisfied for finite size droplets it is sufficient to evaluate n(z, z') for the case of \z\ = \z'\ = r 
. Using the explicit form of the one-particle orbitals Eq. (§) simplifies for this case to 



1 M(N-l) 1 2m 

n(r,r;^) = -e-^ 2 £ (^>_!_ e -^. (11) 

Z7F m=0 Z 

where # is the angle between the z and z' . By making use of the fact that a finite number 
of angular momenta have finite occupation numbers in the Laughlin droplet we obtain an 
exact expression for the AMDF in terms of a sum over a finite number of angles at a single 
radius: 

{m) ~fm(r) M(N-1) + V {lZ) 

where 

1 Jim 

^ - W-^" 1/2 < 13 > 

is proportional to |0 m (r)| 2 and 

Oj = ? 2 "\ ■ (14) 
3 M(N - 1) + 1 v ; 

We evaluate the AMDF numerically by combining a Monte-Carlo evaluation of n(r, r;9j) 
with Eq. (|TJ). 

It is interesting to note that the full AMDF can be determined, and therefore that 
the full one-particle density matrix can be reconstructed, from the angular dependence 
of n(z, z') at any common radius, \z\ = \z'\ = r. This property is unique |TT) to two- 



dimensional systems in the strong magnetic field limit where all electrons are restricted to 
a single Landau level. However, the Monte-Carlo values of n(r, r;9j) will inevitably have 
some statistical uncertainty From Eq. ([12]) we see that the resulting uncertainty in the 
occupation number will be a minimum when r is near the maximum in f m (r) which occurs 
at at r = y/2m = R m . Typical uncertainties in (n m ) become very large unless r is near R m . 
For all the numerical results reported below we estimated (n m ) values from the the angular 
dependence of the density matrix at r near R m . 



The Monte-Carlo calculation evaluates the complex function n(r, r; 9) by a Metropolis (13 
sampling of the positions of particles 2 through N. The weighting factor is 



N 

W(z 2 ,...,z N )= n |^-^r M exp(-^|^| 2 /2). (15) 

Ki<m k>l 



With this weighting factor 

N 



re %e — r,;e lt >j ' 



iV%=l( e - 2 /2JJ {( 

x ( r - rj e^)} M ). (16) 



Qn x />V 



Thus n(r,r;9) is determined up to a constant which is independent of both 9 and r. We 
determine Qn-i/Qn by requiring that the integral of the diagonal elements of the density- 
matrix over position equal N (J dfn(r, r; 9 = 0) = N). n(r, r; 0j) was evaluated in this way 
at a set of radii separated by 0.5£. (n m ) was estimated by averaging results over different 
values of r weighted by the factor e -( r - R m) 2 / 2 . 

Results for the occupation numbers obtained in this way are illustrated in Fig. ([J) for 
M = 3, M = 5, and M = 7 Laughlin droplets with N = 25, N = 20 and N = 15 respectively. 
Monte Carlo calculations were also carried out for smaller size droplets where comparisons 
with analytic results could be made. The two main features apparent in these results are 
the large oscillations in the occupation numbers in the interior of the droplets and the rapid 
decline in the occupation numbers as the outermost orbitals are approached. For the three 
cases illustrated m = 72, 95 and 98 and the occupation numbers for the outermost droplets 
are 1.07 x 10~ 2 , 2.3 x 10 -4 and 2.6 x 10 -5 respectively. In Fig. (|]) we compare the angular 
momentum state occupation numbers near the edge of the droplet with Eq. (|10"D, which we 
find to be accurately satisfied for the droplet sizes studied. We conclude that Luttinger- 
liquid-like behavior should be visible |14| in the low-energy properties even for quite small 



Laughlin droplets. For example for the M = 7, iV = 15, Laughlin droplet m = 98, and the 
Monte Carlo calculation gives (n 93 ) / (n 98 ) = 423 ± 15 compared to the value 462 implied 
by Eq. flT0|). For the droplet sizes illustrated Eq. ( |I~0"D begins to fail badly for k larger than 
~ 10. Away from the edge the occupation numbers reach a maximum value and oscillate 
as the interior is approached. In the limit of extremely large droplets it is known from the 
plasma analogy for Laughlin wave functions that the density deep in the interior approaches 
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(2irM)~ 1 , from which it follows that (n m ) approaches M -1 . Our Monte-Carlo calculations 
show that this limit is approached slowly as N increases. The largest occupation numbers 
occur near the edge of the droplet and the relative excess in this region is larger for larger 
M. For the droplets illustrated the largest occupation numbers are (7150) = 0.52 ± 0.02 for 
M = 3, (n 72 ) = 0.47 ± 0.04 for M = 5 and (n 75 ) = 0.43 ± 0.03 for M = 7. The oscillations 
in occupation number are related to the oscillations in charge density expected near the 



edge of a finite 2D plasma [15j but are much more pronounced since the density averages 
occupation numbers of orbitals with R m near r. 

As we see from the plots of (n m ), the statistical uncertainties, obtained by averaging 
over independent runs, are smaller closer to the edge of the droplet. We believe that this is 
because fewer orbitals contribute importantly to the one-particle density matrix. In fact the 
occupation numbers for orbitals very close to the edge are actually more reliably calculated 
by evaluating the angular dependence of the density-matrix at r substantially larger than 
R mo . In closing we remark that it is in principle possible to determine the occupation num- 
bers uniquely if the radial dependence of the particle-density is known precisely. However, 
this procedure is extremely ill-conditioned and we believe that it is practical only where the 
particle-density is known analytically. We have found it to be impossible to determine the 
AMDF accurately from Monte-Carlo particle densities even for quite small droplets. 

This work was supported by the National Science Foundation under grant DMR-9113911. 
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FIGURES 

FIG. 1. Occupation numbers for M = 3, M = 5, and M = 7, Laughlin droplets with N = 25, 
N = 20, and N = 15 respectively plotted vs. m = The solid line is the particle density in 

(2tt£ 2 )~ 1 units plotted vs.r 2 /2 calculated from the occupation numbers using Eq. (3). The dashed 
line shows the particle density calculated directly; the difference in these two quantities is a measure 
of the error introduced in extracting the occupation numbers from the density matrix. 

FIG. 2. Comparison of Monte Carlo occupation numbers near the edge with the formula 
derived by Wen, Eq. (10). 
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TABLES 

TABLE I. Angular momentum state occupation numbers for some representative small N 
Laughlin droplets 



M, N 


m = 


m = 1 


m = 2 


m = 3 


m = 4 


m = 5 


m = 6 


3,2 


1/4 


3/4 


3/4 


1/4 











3,3 


7/31 


9/31 


18/31 


22/31 


21/31 


12/31 


4/31 


5,2 


1/16 


5/16 


10/16 


10/16 


5/16 


1/16 
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